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ABSTRACT 

We present a sensitive 870 /jm survey of the Extended Chandra Deep Field South (ECDFS) 
combining 310 hours of observing time with the Large Apex BOlometer Camera (LABOCA) 
on the APEX telescope. The LABOCA ECDFS Submillimetre Survey (LESS) covers the full 
30' x 30' field size of the ECDFS and has a uniform noise level of crg7o Mm ~ 1.2 mJy beam -1 . 
LESS is thus the largest contiguous deep submillimetre survey undertaken to date. The noise 
properties of our map show clear evidence that we are beginning to be affected by confusion 
noise. We present a catalog of 126 submillimetre galaxies (SMGs) detected with a significance 
level above 3.7 a, at which level we expect 5 false detections given our map area of 1260 arcmin . 
The ECDFS exhibits a deficit of bright SMGs relative to previously studied blank fields but not 
of normal star- forming galaxies that dominate the extragalactic background light (EBL). This 
is in line with the undcrdcnsities observed for optically defined high redshift source populations 
in the ECDFS (BzKs, DRGs, optically bright AGN and massive K-band selected galaxies). The 
differential source counts in the full field are well described by a power law with a slope of a — 
—3.2, comparable to the results from other fields. We show that the shape of the source counts 
is not uniform across the field. Instead, it steepens in regions with low SMG density. Towards 
the highest overdensities we measure a source-count shape consistent with previous surveys. The 
integrated 870/mi flux densities of our source-count models down to Ssyo/im = 0-5 mJy account 
for > 65% of the estimated EBL from COBE measurements. We have investigated the clustering 
of SMGs in the ECDFS by means of a two-point correlation function and find evidence for strong 
clustering on angular scales < 1' with a significance of 3.4 a . Assuming a power law dependence for 
the correlation function and a typical redshift distribution for the SMGs we derive a characteristic 
angular clustering scale of 0q — 14" ± 7" and a spatial correlation length of = 13 ± 6 Mpc. 

Subject headings: submillimeter: surveys — cosmology: observations — galaxies: evolution — galaxies: 
high-redshift — galaxies: starburst 
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1. Introduction 

One of the most significant findings of the 
IRAS survey was the identification of a popula- 
tion of ultraluminous infrared galaxies (ULIRGs) 
that emit the bulk of their bolometric luminos- 
ity a t far-IR wavelengths ( Sanders fc Mirabel I 
1996). Surveys in the submillimetre and mil- 
limetre wavebands over the past decade have 
shown that ULIRGs are much more common 
at high r edshift compared to the loca l uni- 



result in much more homogeneous coverage and 
mapping, while at the same time allowing for 
a reliable removal of the bright emission from 
the atmosphere. This technique has been used 
at submillime tre and millimetre wavelengths (e.g 
at 35 0/xm, iKovacs et al. 2006; lAustermann et al 



2009, at llOO^m), however, no deep survey have 



verse Ce.g. iBareer etal.1 Il999t 
Borvs et all I2003L IWebb et al ' 



.. .... .. .Greye et al 

2004t lLaurent et all 120051: ICoppin et all 12 006: 



Cowi e et al 
200.4 



2002t 



Pope et alj|2006t iBertoldi et al.ll2007t iBeelen et aT 



2008tlKnudsen et al.l2008HScott et aU2008HAust 



2009) . These surveys show that the comoving vol- 
ume density of luminous submillimetre galaxies 
(SMGs) increases by a factor of 1000 out to z ~ 2 
(jChapman et al.l 120051 ) . Therefore luminous ob- 
scured galaxies at high redshift could dominate 
the tot al bolometric emission of galaxies at th ose 
epochs (|Blain et al.lll999t iLe Floch etaDl2005l ) 

The identification and study of submillimetre 
galaxies has proved challenging since their first 
detection. The limited mapping speed of typi- 
cal (sub) millimetre bolometer cameras meant that 
only few very bright examples have been found, al- 
though gr avitational lensing initially aided some- 
what (e.g. lSmail. Ivison fc Blamlll997t Ivison et al 



1998). Attempts to map large fields at submillime- 
tre wavelengths have involved the use of patch- 
works of small "jiggle" maps (e.g. ICoppin et al 



2006) or mixtures of single-bolometer photom- 
etry, small jiggle maps and shallow scan maps 
used to construct a "Super- map" of GOODS-N 
(|Borvs etal.ll2003t|Pope et~ai1l2006f) . Both of these 
approaches raise concerns about the homogene- 
ity of the resulting maps and hence the relia- 
bility of the resulting source catalogues. Scan 
maps, where the array is continuously moved on 
the sky to trace out a closed pattern, should 
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employed such a technique in the 870-/^m window 
where most of the published work on SMGs has 
been undertaken. Drawing this distinction be- 
tween 870-/im and 1100-/um surveys may appear 
surprising given the modest difference between the 
two wavelengths and the assumed unstructured 
nature of the dust spectrum at these wavelengths. 

Despite only a 25% difference in the two wave- 

' rraaml e ^l§rb^ ths, there are hints of significant differences 
in the po pulations identified at 870- urn and 1100 - 
fj,m (e.g. iGreve et al.l 120041 : lYounger et all l2008h . 
although these may in turn reflect the different 
mapping techniques used in individual studies. 

The advent of the n ew Large APEX B olome- 
ter Camera (LABOCA. ISiringo et aT1l2009l ). with 
an instantaneous 11 . 4' field of view, o n the 12- 
m APEX telescope ( Glisten et al.l 12006^ provided 
the opportunity to undertake the first sensitive 
and uniform panoramic survey of the extragalactic 
sky at 870/^m. To exploit this opportunity a num- 
ber of groups within the Max Planck Gesellschaft 
(MPG) and the European Southern Observatory 
(ESO) communities proposed a joint public legacy 
survey of the Extended Chandra Deep Field South 
(hereafter ECDFS) to the MPG and ESO time 
allocation committees: the LABOCA ECDFS 
submillimetre survey (LESS). The ECDFS cov- 
ers a 0.5° x 0.5° region centered on the Chandra 
Deep Field South (CDFS) at RA 03 /l 32 m 28 s .O 
Dec -27°48'30".0. This field has very low far- 
infrared backgrounds and good ALMA visibility 
and hence has become one of the pre-eminent 
fields for cosmological survey science. As a result, 
the ECDFS is unique in the Southern Hemisphere 
in the combination of area, depth and spatial res- 
olution of its multiwavelength coverage from X- 
rays through optical, near- and mid-infrared to 
the far-infrared and radio regimes. The central 
part of this field is coincident with the CDFS 
( Giacconi et al. 2002f) which has now reached a 
depth of 2 Ms (|Luo et all l2008f) and the deep 
Hubble Space Te lescope (HST) i magin g of the 
GOODS-S field (|Giavalisco et all l2004h and the 
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Hubble Ultra Deep Field (UDF, iBeckwith et al 
20061 ). In addition to the extremely deep observa- 
tions of the central regions of this field as part of 
the CDFS, GOODS and Hubble Ultra Deep Field 
surveys, the full 0.5° field has extensive multiwave- 
length imaging available including: 250-ks Chan- 



dra in tegrations over the whole field (jLehmer et al.l 
20051 ): deep and multi-band optical imaging by 
COMBO- 17 dWolf et alJ l20oi 2008) and MUSYC 
(|Gawiser et al.l 120061) ncluding HST imag ing for 
the GEMS project (ICaldwell et al.l 120081): near - 
infrared imaging by MUSYC (|Tavlor et alj|2009h : 
deep mid- infrared imaging w ith IRAC as part of 
SIMPLE (|Damen et alj|2009h and using the MIPS 
instrument at 24, 70 and 160/xm by FIDEL (Dick- 
inson ct al. in prep.). Longer wavelength coverage 
comes from BLAST (jPevlin et al.l 120091) at 250, 
350 and 500/im (and in the near future from Her- 
sc hel) , while radio cov erag e of this field i s repo rted 
bv lMiller etal] (|2008l ) and llvison etail (|2009h . 

The LABOCA survey of the ECDFS adds a 
waveband that pin-points the thermal emission 
from luminous dusty galaxies at z ~ 1-8: a power- 
ful addition to this singularly well-studied region - 
ideally placed for VLT observations and early sci- 
ence follow-up with ALMA. The completed LESS 
project provides a representative, homogeneous 
and statistically-reliable sample of the SMGs with 
the high-quality, multiwavelength data required 
to yield identifications, constrain their redshifts, 
bolometric luminosities and power sources and 
hence determine their contribution to the total 
star formation density at high redshift. These 
sources can be related in unprecedented precision 
to other populations of AGN and galaxies within 
the same volume to understand the place of SMGs 
in the formation and evolution of massive galax- 
ies at high redshift. The survey is also sufficiently 
large that it should also yield examples of rare 
classes of SMGs, such as ve ry high redshift sources, 



z > 4 (|Coppin etal.ll2009l ). These same data also 



provide submillimetre coverage of large numbers 
of high-redshift galaxies and AGN to determine 
their bulk submillimetre properties from the stack- 
ing analysis of sub-samples as a func tion of pop- 



ulation, redshift, environment, etc. ( Greve et al 



20091: iLutz etal.ll2009l) . Together, these two tech- 
niques allow us to sample two orders of magnitude 
in bolometric luminosity - from hyperluminous in- 
frared galaxies with 10 13 L0 which are directly de- 



tected in the LABOCA maps, down to luminous 
infrared galaxies at 10 11 L Q which are detected 
statistically through stacking. This range in lu- 
minosity encompasses the variety of populations 
expected to dominate the bolometric emission at 
z ~ 1-3 and the cosmic submillimetre background. 

In this paper we present a detailed description 
of the observations, reduction and analysis of the 
LABOCA observations of the ECDFS and the re- 
sulting catalogue of submillimetre galaxies. The 
observations are described in §2, §3 presents our 
results and we discuss these in §4. Finally, in 
§5, we give our summary and the main conclu- 
sions of this work. We assume a cosmology with 
H = 70kms- x Mpc" 1 , fl A = 0.7 and Q M = 0.3. 

2. Observations and data reduction 

Observations were carried out usi ng the Large 
APE X Bolometer Camera LA BOCA dSirineo et al. 
20091 ) on the APEX telescope (|Giisten et al .1120061 ) 
at Llano de Chajnantor in Chile. LABOCA is an 
array of 295 composite bolometers with neutron- 
transmutation-doped (NTD) germanium thermis- 
tors. The bolometers are AC-biased and operated 
in total power mode. Real-time signal processing 
of the 1 kHz data stream includes digital anti- 
alias filtering and down-sampling to 25 Hz. The 
radiation is coupled onto the detectors through 
an array of conical feed horns whose layout leads 
to a double beam spaced distribution of the in- 
dividual beams in a hexagonal configuration over 
the 11'. 4 field of view. The center wavelength of 
LABOCA is 870^im (345 GHz) and its passband 
has a FWHM of - 150/^m (60 GHz). The mea- 
sured angular resolution of each beam is 19". 2 
FWHM. 

The observations^ were carried out between 
May 2007 and November 2008 in mostly excellent 
weather conditions with an average precipitable 
water vapor (PWV) of 0.5 mm corresponding to a 
zenith opacity of 0.2 at the observing wavelength. 
The mapping pattern was chosen to give a uni- 
form coverage across a 30' x 30' area centered at 
RA 03 /l 32 m 29 s .O Dec -27°48'47".0. Mapping was 
performed by alternating rectangular, horizontal 
on-the-fly (OTF) scans with a raster of spirals pat- 
tern. OTF maps were done with a scanning veloc- 



iProgramme IDs 078.F-9028(A), 079.F-9500(A), 080.A- 
3023(A) and 081.F-9500(A). 
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ity of 2 arcmins -1 and a spacing orthogonal to the 
scanning direction of 1'. For the spiral mode, the 
telescope traces in two scans spirals with radii be- 
tween 2' and 3' at 16 and 9 positions (the raster) 
sp aced by 10 / in azimu th and elevation (see Fig. 9 



Siringo et al. (l2009h for a plot of this scanning 



pattern). The radii and spacings of the spirals 
were optimized for uniform noise coverage across 
the 30' x 30' region, while keeping telescope over- 
heads at a minimum. The scanning speed varies 
between 2 - 3arcmins _1 , modulating the source 
signals into the useful post-detection frequency 
band (0.1 to 12.5 Hz) of LABOCA, while provid- 
ing at least 3 measurements per beam at the data 
rate of 25 samples per second even at the highest 
scanning velocity. 

Absolute flux calibration was achieved through 
observations of Mars, Uranus and Neptune as well 
as secondary calibrators (V883 Ori, NGC 2071 and 
VYCMa) and was found to be accurate within 
8.5% (rms). The atmospheric attenuation was de- 
termined via skydips every ~ 2 hours as well as 
from independent data from the APEX radiome- 
ter which measures the li ne of sight water va por 
column every minute (see ISiringo et al.l 20091 . for 
a more detailed description). Focus settings were 
typically determined a few times per night and 
checked during sunrise depending on the availabil- 
ity of suitable sources. Pointing was checked on 
the nearby quasars PMNJ0457-2324, PMNJ0106- 
4034 and PMNJ0403-3605 and found to be stable 
within 3" (rms). 

The data were reduced using the Bolometer ar- 
ray data Analysis software (BoA, Schuller et al. in 
prep.). Reduction steps on the time series (time 
ordered data of each bolometer) include tempera- 
ture drift correction based on two " blind" bolome- 
ters (whose horns have been sealed to block the 
sky signal), flat fielding, calibration, opacity cor- 
rection, flagging of unsuitable data (bad bolome- 
ters and/or data taken outside reasonable tele- 
scope scanning velocity and acceleration limits) 
as well as de-spiking. The correlated noise re- 
moval was performed using the median signal of 
all bolometers in the array as well as on groups 
of bolometers related by the wir ing and in the 
electronics (see Siringo et al. 20091 ). After the de- 
correlation, frequencies below 0.5 Hz were filtered 
using a noise whitening algorithm. Dead or noisy 
bolometers were identified based on the noise level 



of the reduced time series for each detector. The 
number of useful bolometers is typically ~ 250. 
The data quality of each scan was evaluated using 
the mean rms of all useful detectors before correct- 
ing for the atmospheric attenuation (which effec- 
tively measures the instrumental noise equivalent 
flux density, NEFD) and based on the number of 
spikes (measuring interferences). After omitting 
bad data we are left with an on-source integra- 
tion time of ~ 200 hrs. Each good scan was then 
gridded into a spatial intensity and weighting map 
with a pixel size of 6" x 6". This pixel size (~l/3 
of the beam size) well oversamples the beam and 
therefore accurately preserves the spatial informa- 
tion in the map. Weights are calculated based on 
the rms of each time series contributing to a cer- 
tain grid point in the map. Individual maps were 
coadded noise-weighted. The resulting map was 
used in a second iteration of the reduction to flag 
those parts of the time streams with sources of a 
signal to noise ratio > 3.7o\ This cut off is de- 
fined by our source extraction algorithm. The re- 
duction with the significant sources flagged guar- 
antees that the source fluxes are not affected by 
filtering and baseline subtraction and essentially 
corresponds to the very same reduction steps that 
have been performed on the calibrators. 

To remove remaining low frequency noise arte- 
facts we convolved the final coadded map with a 
90" Gaussian kernel and subtracted the resulting 
large scale structures from the unsmoothed map. 
The convolution kernel has been adjusted to match 
the low frequency excess in the map. This step is 
effectively equivalent to the low frequency behav- 
ior of an optimal point-source (Wiener) filtering 
operation (jLaurent et al.l 120051 ). The effective de- 
crease of the source fluxes (~ 5%) for this well 
defined operation has been taken into account by 
scaling the fluxes accordingly. Finally the map 
was beam smoothed (convolved by the beam size 
of 19.2") to optimally filter the high frequencies for 
point sources. This step reduces the spatial res- 
olution to w 27". The signal and signal to noise 
(S/N) presentations of our final data product is 
shown in Fig.[TJ 

To ensure that aboves reduction steps do not 
affect the flux calibration of our map we per- 
formed the same reduction steps on simulated time 
streams with known source fluxes and artificial 
correlated and Gaussian noise. These tests veri- 
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Fig. 1. — Flux (top) and signal to noise (bottom) map of the ECDFS at a spatial resolution of 27" (beam 
smoothed). The white box shows the full 30' x 30' of the ECDFS as defined by the GEMS project. The 
white contour shows the 1.6 mJy beam -1 noise level that has been used to define the field size for source 
extraction yielding a search area of 1260 sq. arcmins. The circles in the top panel indicate the location of 
the sources listed in Tableful 
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Fig. 2. — Top: Normalized flux map of 
PMNJ0403-3605 (flux 590 mJy), the most fre- 
quently used pointing source near the ECDFS, re- 
duced in the same way as the ECDFS scans. Con- 
tours are shown at 5, 15 (black) and 30, 50, 70 and 
90% (white) of the peak flux density. Bottom: Ra- 
dial averaged beam profile. The solid line shows a 
Gaussian fit that yields a FWHM of 19". 2. 

fled that our calibration scheme is accurate 
to ~ 5%. Furthermore we reduced data of 
PSS 2322+1944, a z = 4.1 QSO which has 
been observed during the science verification 
of LABOCA, in the same way as the ECDFS 
data. From this measurement we find SWo/otj = 
21.1±2.5mJy in good agreeme nt with the SCUB A 
measurements (22.5 ± 2.5 mJy. Ilsaak et ah 2002f h 
Finally we determined the point spread func- 
tion (PSF) in our map by applying the same data 
reduction steps to PMNJ0403-3605, the most fre- 
quently used pointing source near the ECDFS. 
The beam profile (before beam smoothing) is 



shown in Fig.[2]and is well described by a Gaussian 
with FWHM of 19.2". The faint negative struc- 
ture at radii between 25" and 45" is due to the 
combined effect of the correlated noise removal, 
the low-frequency filtering on the time series and 
the spatial large scale filtering. 

3. Results 

3.1. Noise Properties 

To investigate the noise properties of our 
LABOCA map we have created 100 pure noise re- 
alizations of the data by randomly inverting half of 
the maps of individual scans during the coadding 
(e.g. IPerera et al.ll2008t fecott et al.ll2008f) . All im- 
age processing steps match the real map including 
the large scale filtering. These so-called "jack- 
knife maps" are therefore free of any astronomical 
signal and at the same time represent the noise 
structure of the data. 

We first investigate how the noise in our map 
integrates down with time. This is shown in Fig. [3] 
where we plot the rms noise level measured on the 
central 30 x 30 of our map as a function of the 
integration time per beam. We have generated 
down-integration curves for the pure noise realiza- 
tions, the unmodified data and for the data after 
subtracting the full source catalog (see Sect. 13.2.4]) 
from each scan. All three computations were done 
at the original spatial resolution as well as for 
the beam-smoothed data. Ten computations with 
randomized scan order were performed for each 
method. Fig. [3] shows the average of these compu- 
tations. 

As expected, the pure noise realizations inte- 
grate down proportionally to y/TJt independent 
of the spatial resolution. At the original spa- 
tial resolution of 19.2" the influence of the source 
signals is small and even the down-integrating of 
the data including all sources follows this behav- 
ior very closely. Some excess noise is barely vis- 
ible after several hundred seconds of integration 
time per beam. The noise excess becomes very 
pronounced for the beam smoothed data and it 
remains visible even after the subtraction of the 
source catalog from the data. In fact, much of the 
noise access (65%) remains in our beam-smoothed 
down-integrating curves after source subtraction. 
This shows that the noise in our map is limited by 
sources fainter than our catalog limit (~ 4.7 mJy, 
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Noise (Jackknifes) 

Data w. sources subtracted 

Data 




10 2 10 3 
integration time/beam [s] 

Fig. 3. — Average noise level in the ECDFS as 
a function of integration time. The upper lines 
represent the unsmoothed data at a spatial res- 
olution of 19.2". The three lines show the rms 
noise vs. integration time for the pure noise maps 
(jackknifes - solid line), the flux map after subtrac- 
tion of all sources in the catalog (see Sect. l3.2."4l - 
dashed line) and for the flux map including all 
sources (dotted line). The lower curves are the 
same but for beam-smoothed data (27.2" spatial 
resolution). The grey solid line has a slope of -0.5 
and shows the expected behavior for pure noise. 

Sect. l3.2.^|) or in other words that our map starts 
to be confusion limited. We have estimated the 
confusion noise arising from these faint sources by 
fitting the down integration curve for the beam 
smoothed data using 

a ohs (t) = y/alit) + al (1) 

where <7 b s , o~ n an d cr c are the observed, instrumen- 
tal/atmospheric and confusion noise terms respec- 
tively. From this fit we derive a confusion noise of 
cr c w 0.9 mJy beam" 1 at 27" resolution. We note 
that this level of confusion noise is consistent with 
the simulated down-integration curves based on 
our source counts (see Sect. l3.3.Tj) . 

Given that our map rms is affected by confu- 
sion noise, the question arises which noise level is 
appropriate for computing the noise distribution 
across the field and the corresponding signal-to- 




-5 5 10 15 

Pixel flux [mjy/beam] 




Pixel S/N 

Fig. 4. — Top: Observed beam smoothed flux his- 
togram (black) compared to the average flux his- 
togram of the 100 jackknifes (normalized to the 
same peak) in the ECDFS. The positive access is 
due to sources, the broadening of the negative part 
of the observed flux histogram is due to the confu- 
sion noise (see text). Bottom: Beam smoothed 
signal-to-noise (S/N) histogram computed from 
the scaled weights of the map resulting in a Gaus- 
sian with a — 1 (see Sect. 13. ip . 

noise (S/N) map, which is used for source extrac- 
tion. Usually, the weights of the data (reflecting 
the rms weighted integration time in each pixel) 
are used to derive the noise distribution. This ap- 
proach, however, is equivalent to using the jack- 
knife noise, neglects the confusion noise and there- 
fore overestimates the S/N ratio. We here use 
the noise based on the weights with a scaling to 
take the confusion noise into account. The scaling 
was determined from the flux and S /N histograms 
in the ECDFS. The basic principle is shown in 
Fig. S] (top) where we compare the average beam 
smoothed flux pixel histograms of the jackknifes 
to a histogram of the real map. The figure shows 
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Jy/b 



Fig. 5.— Noise map of the ECDFS at 27" resolu- 
tion (beam smoothed). The circles show the loca- 
tion of the 126 sources listed in TableO The black 
box shows the full 30' x 30' of the ECDFS. The 
white contour shows the field size that was used 
for source extraction. The black contours show the 
noise level at 1.1, 1.2, 1.4 and 1.6 mJybeam -1 . 



that the negative part of the observed histogram is 
significantly broader than the jackknife histogram. 
This is because inserting sources in a pure noise 
map (e.g. the jackknife) will not only result in 
the positive flux tail, but also shift the entire 
flux histogram to positive values and broaden the 
Gaussian part of the distribution (because it is no 
longer centered around zero and sources also fall 
on the formally negative part of the flux distribu- 
tion). As the zero point (true total power infor- 
mation) is undetermined in our reduction (e.g. via 
baseline subtraction) , the observed flux histogram 
is roughly centered at zero and only the broaden- 
ing remains an observable compared to the jack- 
knife. This broadening can be used to take the 
confusion noise into account. In practice, we have 
scaled the rms map derived from the weights such 
that the negative part of the Gaussian signal-to- 
noise histogram (Fig. 2] bottom) has a a of unity. 

The resulting spatial distribution of the beam- 
smoothed noise level across the imaged field is 
shown in Fig. [5] The deepest part of our map 
has a noise level of 1.10 mJybeam -1 , the typi- 
cal rms at the edges of the 30 x 30 field is 1.25 
mJybeam -1 . The average noise level across the 
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Fig. 6. — Area and mean 850/im noise level for se- 
lected mm/submm deep fields compared to LESS. 
Filled and open symbols represent 850/mi and 
1.1mm surveys respectively. For the later we 
have scaled the noise level by a factor of 3.8, the 
mean flux ratio of SMGs between both wavelength 
(|Greve et al.ll2008[h 



field is 1.17 mJybeam -1 with a dispersion of 40 
/iJybeam -1 only. For comparison with previous 
work we compare the area and the noise level of 
selected mm/submm deep field survey to the LESS 
in Fig.[6j 

3.2. Source Catalog 

3.2.1. Source extraction algorithm 

From the beam-smoothed map we have ex- 
tracted sources using the false detection rate al- 
gori thm (for a descrip tion of this method see 
e.g. iHopkins et al.l l2002l) of the CRUSH package 
(iKovacsfeOoiT The choice of beam smoothing be- 
comes obvious if one considers the maximum likeli- 
hood amplitude A of a point source (i.e. the beam 
B) fitted at a given position x in the map <S(x). 
Consider a weighted \ 2 of the fit defined in the 
usual way as: 

X 2 = ]T w(x') [S(*) A(x)B(x - x')] 2 . (2) 



S 



Then, the \ 2 minimizing condition d\ 2 /dA = 
yields a maximum-likelihood amplitude: 



in Tabled 



Position uncertainly 



A(x) 



E x ,m(x')2?(x-x')S(x') 
J] x ,«;(x')S 2 (x-x') 



(3) 



which is effectively the weighted beam smoothed 
image. Thus, the flux values of the beam- 
smoothed map essentially measure the fitted 
beam amplitudes at each map position (see also 
Serjeant et alll2003t iGawiser et al.ll2006l ). 

The algorithm uses as in input parameter the 
allowed number of false detections. Based on this 
number and the field/beam size we calculate as 
an initial search criterion a detection signal-to- 
noise cutoff level assuming Gaussian noise statis- 
tics. When identifying source candidates, we al- 
low for the possibility that the true peak may fall 
between pixels by appropriately relaxing the ini- 
tial search signal-to-noise ratio. Then for each 
source candidate, we interpolate the neighboring 
pixel values to estimate the underlying peak and 
its position, and keep only candidates that meet 
the original detection S/N level. We also apply the 
same procedure to identify negative noise peaks. 
These we use to check if the false detection rate 
is consistent with Gaussian noise, or to adjust the 
empirical noise distribution for possible deviations 
thereof. 

The source candidates thus identified are re- 
moved from the map using the appropriate large- 
scale-structure (LSS) filtered beam profiles. LSS 
filtering corrections are also calculated for each 
pair of source candidates. The maps are flagged 
around the extracted source positions. After the 
extraction step, the zero flux level of the map is 
re-estimated via the mode of the pixel distribu- 
tion (which we consider the most robust measure 
in this case). We also re-estimate the width of 
the noise distribution using robust measures. The 
extraction steps are repeated until no further can- 
didates are identified. 

For each candidate, we calculate a detection 
probability based on the number of sources de- 
tected beforehand and of the number noise peaks 
found (or expected) below the corresponding in- 
verted significance level, under the assumption of a 
symmetric noise probability distribution. Finally, 
the candidates are sorted in order of decreasing 
detection probability. For each source we indicate 
the corresponding cumulative false detection rate 
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Fig. 7. — Top: Position uncertainty as a func- 
tion of observed signal-to-noise ratio. Error bars 
show the rms of the extracted positions in each 
signal-to-noise bin, the red line shows the ex- 
pected posit i onal u ncertainties using Eq. B22 from 



Ivison et alj (|2007| ). Middle and bottom: Com- 
pleteness and flux boosting as a function of intrin- 
sic source flux determined from Monte-Carlo simu- 
lation. The dashed vertical line in the lower panel 
shows the lowest deboosted flux in our source cat- 
alog. 



3.2.2. Completeness and Position accuracy 

To test the reliability of our extraction process 
and to obtain information on the completeness of 
our catalog we ran our extraction algorithm on 
100 jackknife noise realizations after inserting ar- 
tificial sources. The sources were added at random 
positions for each jackknife map and therefore rep- 
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resent an unclustered source population. For the 
source flux distribution we used a power law con- 
sistent with the differential source counts derived 
in Sect. 13. 31 Simulated sources were inserted down 
to a flux level of 1.0 mJy, i.e. well below our de- 
tection threshold for individual source. The total 
number of simulated source per jackknife map was 
~ 6000. The large number of sources implies that 
the effect of confusion noise is taken into account 
in our extraction. The completeness computed 
from these simulations is shown in middle panel 
of Fig. [3 The figure shows that our source extrac- 
tion is complete (> 95%) down to an input source 
flux of ~ 6.5 mJy while the 50% completeness level 
occurs at ~ 4.0 mJy. 

We also used our Monte-Carlo simulations to 
estimate the accuracy of the source coordinates 
determined by our source extraction algorithm. 
This is shown in the top panel of Fig. [7J where we 
plot the positional difference between input and 
extracted position as a function of the measured 
signal-to-noise ratio of the extracted sources. The 
mean positional un certainty is i n good agreement 
with Eq.B22 from llvison et al.1 (|2007l red line). 
The scatter in the low S/N bins (< 5), however, 
is large which implies that the positional uncer- 
tainties from aboves equation are only correct in 
a statistical sense and that the true offset of an 
individual source can be much larger. From the 
simulations we find that for our extraction limit 
(S/N > 3.7) > 95% of the sources will have a po- 
sitional accuracy better than 8". For the brightest 
sources (S/N > 7) the 95% confidence radius is 
» 3.5". 

We note that because of the large number of in- 
dependent observations in our map the telescope 
pointing accuracy will not affect the astrometry 
of our sources but result in a small (< 3%) spa- 
tial smearing of their signals. The accuracy of the 
absolute astrometry of our map has been tested 
based on our stacking analysis of a sampl e of near- 



IR selected galaxies (jGreve et al. 2009). These 



stacks give signals with up to 20er significance cen- 
tered on the expected position for different K-band 
selected source populations and result in signal 
profiles in agreement with the beam shape. From 
this we conclude that there is no evidence for an 
overall mean astrometric error of our data. This 
is also confirmed by a comparison of the 1.4 GHz 
radio relative to the submm positions (Ivison et al. 



in prep.). 

3.2.3. Flux boosting correction 

Signal-to-noise limited source catalogs carry a 
selection bias from an overabundance of sources 
whose ap parent flux is positively enhanced b y 
noise (e.g. lHogg fc TurneJll998llScott et al.ll2002h . 
This is shown in the bottom panel of Fig. [7] where 
we compare the extracted flux densities to the 
input flux densities of our Monte-Carlo simula- 
tions. The discrepancy between intrinsic fluxes 
and their detection values becomes noticeable be- 
low ~ 6 mJy. 

This effect arises because the steepness of 
submm number counts implies that an observed 
flux Sobs more often arises from intrinsic fluxes 
S < Sobs and less often from S > Sobs- The av- 
erage flux value (S) behind an observed flux S oos 
can be calculated statistically: 



(S) 



S p{S\S ohs ) dS 



(4) 



According to Bayesian theory, p(Si\S ha) 
p{S bs\S) p{S). The probability density p(S bs\S) 
of observing a flux S oos for an underlying flux S is 
simply the noise distribution evaluated at the flux 
difference, i.e. n(S — S oos ). For Gaussian noise 
n(x) is calculated as (aV^n) -1 exp(— x 2 /2a 2 ). 
The term p{S) is the probability density of flux S 
in the underlying noiseless flux distribution of the 
map, which accounts for the possibility of over- 
lapping sources. As such, p(S) is given by Eq. [7] 
and is a direct product of our P(D) analysis for 
the source number counts (see Sec. 13.3.11) . 

A slight complication arises because the map 
zeroing is biased by the the presence of sources 
below detection level. If the map zero level corre- 
sponds to an intrinsic flux SS, then a map flux S 
really belongs to an underlying value S+SS. For- 
tunately, this SS is also readily produced by the 
P(D) analysis. With the necessary modifications 
in place we can calculate deboosted fluxes using 



(S) =ss + 



J S n(S — S bs) p(S + SS) dS 
J n(S — Sobs) p(S + SS) dS 



(5) 



The integration can be performed numerically. 
The corresponding uncertainty a of the deboosted 

\s-(s)f), 



flux can also be calculated as a 2 
i.e., as 
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2 _ f{S + 6S- (S)f n(S - Sobs) p(S + SS) dS 
Jn(S-S ohs )p(S + SS)dS 

(6) 

The uncertainty of the deboosted flux is typ- 
ically larger than the measurement uncertainty. 
The effect is more more pronounced at the lower 
fluxes. 

To calculate the deboosted fluxes for each 
source we have used p(S) and SS from the sin- 
gle powerlaw fit from our P(D) analysis derived 
in Sect. 13. 3. Tl Note that this is not an iterative 
process because the P(D) fitting does not require 
information on the underlying counts. 

3.2.4- The LESS source catalog 

While our map has a uniform noise of 1.2 mJy 
beam -1 over 900 sq. arcmin, an additional 360 
sq. arcmin has only slightly higher noise, < 1.6 
mJybeam -1 (i.e. better than most previous sur- 
veys, see Fig.[6j, hence we expand our search area 
to allow us to find slightly brighter sources out- 
side the uniform region (see Fig. [I]). To construct 
a robust catalog we restricted the extraction in 
such a way to give statistically 5 false detections. 
This yields a total of 126 sources that are listed in 
Table H The false detection rate (FDR) of 5 im- 
plies a signal to noise cut above 3.7a. The faintest 
sources in our catalog have measured flux densi- 
ties of 4.6 mJy, at which level the completeness is 
> 70% (Fig. EJ|. In the Table we order the in- 
dividual sources by their signal-to-noise ratio and 
we give their IAU name, the source position, the 
measured source flux with the map noise as uncer- 
tainty, the de-boosted source flux and its uncer- 
tainty using a Bayesian approach (see Sect. 13.273)) , 
the signal-to-noise ratio and the expected number 
of false detections for all sources including the cor- 
responding entry in the Table. From the last entry 
in the Table it can be seen that only 10 additional 
sources are included in the catalog if we increase 
the FDR from 3 to 5 for the full catalog. This im- 
plies that 20% of the additional sources are likely 
to be false. Therefore deeper source extractions 
would not yield reliable information. 

3.2.5. Tests on the LESS source catalog 

To test the reliability of our source catalog we 
have compared the FDR from our Monte-Carlo 



simulations to the FDR expectation for our source 
extraction algorithm. A source was considered to 
be detected if the extraction fell within a 8" search 
radius (the maximum positional uncertainty ex- 
pected for our extraction, see l3.2.2[) from the input 
position. The FDRs derived in both ways agree 
very well with a slight tendency of our extraction 
algorithm to overestimate the FDR (for 5 expected 
false detection we find 4 from our Monte-Carlo 
simulations). 

To verify that our map does not contain false 
sources due to artefacts in the data or the data 
processing we split our observations into two 
parts. This was done by splitting the randomized 
scan list (2370 scans in total) into two lists with 
roughly equal integration time yields two indepen- 
dent maps with noise levels of « 1.7 mJy beam -1 . 
We then performed our source extraction on both 
maps using the same significance level we used 
for the generation of our source catalog (3.7er). 
This yields 59 and 60 extracted sources for map 
1 and map 2, respectively which is in excellent 
agreement with the expected numbers considering 
the y/2 increase of the noise level and the number 
counts derived from the full map. 

From these extractions we find that 22 sources 
in the LESS catalog are detected in both submaps 
and that ~ 70% of the sources extracted from both 
submaps are also in the LESS catalog. In the last 
entry of Tabled] we have indicated for each source 
of the LESS catalog whether they are detected in 
one, both or none of the individual maps. Note, 
that those sources extracted from both submaps 
which are not in the LESS catalog are not neces- 
sarily false detections as the different noise struc- 
ture of the submap may boost other faint sources 
above the extraction threshold. 

A comparison of the sources detected on map 
1 and 2 yields an overlap of both catalogs of 
~ 40%. For those sources that only appear in 
one of the two catalogs (37/38 for map 1/2) we 
have extracted the SNR peaks in a 8" search ra- 
dius on the map where the source is not detected 
and computed the probability that no source is 
present in this aperture. From this analysis we 
find that the probability of a false detection ex- 
ceeds 10% for just 7(/6) source for submap l(/2) 
which is in reasonable agreement with the 5 false 
detections expected from our extractions. We fur- 
thermore analyzed the stacked signal at the po- 
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Pixel flux [mJy/beam] 

Fig. 8. — Top: Histogram of the flux density dis- 
tribution of all pixels in the central 30' x 30' af- 
ter subtracting all sources from our catalog (Table 
[2|) . The curve represents a Gaussian fit to the un- 
derlying noise histogram fitted between -5.0 and 
1.0 mJy. The region shows the residual positive 
access. Bottom: Linear presentation of the flux 
excess. Errors are 3 a Poisson uncertainties. 



sitions of the non-detec tions. For t his w e used 
the method described bv lGreve et all tl2009h which 
yields « 12 a detections for both submaps with an 
average flux density of 3.5 ± 0.3 mJy for all posi- 
tions and 4.3 ± 0.35 mJy if we exclude the 7(/6) 
false detections. 

Finally we have investigated the stacked sig- 
nal in maps 1 and 2 at each position at which a 
source is extracted in the LESS catalog, but not in 
both submaps. Again the stacked signals in both 
submaps give very similar results with ss 20 a de- 
tections and an average flux of 4.2 ± 0.2 mJy. The 
intensity of stacked signals is in good agreement 
with the mean deboosted source flux of the LESS 
sources entering into the stack for both maps. 
These tests show, that our FDR extraction yields 
reasonable results and that the average fluxes of 
the faint sources in the LESS catalog can be repro- 
duced in the stacking signals of two independent 
submaps. 



3.3. Differential source counts 

Most studies that have addressed mm/submm 
number counts from deep-field surveys have used 
the extracted sources, e.g. the resulting source 
catalogs to fit models to the differential or in- 
tegrated nu mber counts as a function of flux 
density (e.g. iBarger et al Il999l: iBlain et aD 19991 
Scott et all 120021: iGreve et all 120041: ICoppm et al 



20061: iKnudsen et all 20081 : Perera et alj 1200a 
Austermann et al.l 12009 ). This approach relies on 
the correct determination of the completeness as 
well as on the flux de-boosting in order to extract 
the required information. In particular the latter 
step is problematic because the de-boosting it- 
self requires information on t he underlying sourc e 



count distribution (see e.g. Coppin et al. 20051) . 



Furthermore, blank field surveys contain also in- 
formation on sources fainter than th e typical cutoff 



levels used to extract sources (e.g. iPeacock et al 
2000). This is shown for our LABOCA data in 



Fig. [8] where we show the flux pixel histogram 
after subtracting all sources in our source cata- 
log from the map. The diagram shows a signif- 
icant access of pixels with positive flux densities 
which is due to sources below our source extrac- 
tion limit. The increasing number of noise peaks 
at the same flux density level prevents an ex- 
traction directly from the map without increasing 
significantly the number of false detections in the 
analysis. One may, however, ignore the position 
information completely and derive the underlying 
source count distribution directly from the flux 
pixel histogram through a P(D) analysis (see e.g 
Condon I Il97i iHughes etai1ll998t iMalonev etal 
20051) . In the bottom panel of Fig. [5] we show the 



pixel histogram of the positive excess after sub- 
tracting a Gaussian noise distribution. The latter 
has been determined by a fit to the negative part 
of the flux histogram. The errors in the figure 
are 3 a Poisson uncertainties. From the figure it 
can be seen that the excess becomes significant for 
flux densities above ~ 1.5 mJy which implies that 
the source count distribution can be extracted to 
a much fainter level than the limits inherent to 
the source extraction method. 

3.3.1. P{D) analysis 

Our approach is to directly use the informa- 
tion from the observed flux distribution in our 
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Crete bins i, each containing Ni sources with flux 
Si, then the characteristic function it (i.e. Fourier 
transform) of the compound distribution is: 



Fig. 9. — Top: Logarithmic presentation of the 
observed flux density histogram (black) and fit re- 
sults from the P(D) analysis (red) in the central 
30' x 30' of the E-CDFS. Bottom: Linear presen- 
tation of the residual between the observed and 
fitted flux density histogram. Errors are 3 er Pois- 
son uncertainties. The grey histogram in the top 
panel shows for comparison the flux density distri- 
bution for a Schechter function as derived for the 
SHADES survey with N' matched to the LESS 
source counts (dashed grey line in Fig. [TOl left). 

map. Were there just a single source with flux 
S present at some random location in a noiseless 
map, it would give rise to a probability distribu- 
tion P(y\S) of pixel values y. If the map contains 
several such randomly distributed sources with 
fluxes Si...S n , then the probability P(y\Si...S n ) 
of observing a value y (in flux or S/N) can be 
related to the individual probability distributions 
P(y\Si) arising from each individual source and 
the measurement noise distribution n(y) as, 

P(y\S!...S n ) = P(y\S 1 )®...®P(y\S n )®n(y), (7) 

which simply follows the compounding rule for in- 
dependent random variables via convolution. As 
such, it can be rewritten as the product of the 
complex Fourier transforms (P f± it, and of n f± 
v) of the distributions: ~ir(S\...S n ) = tt{S\) ■ 
ir(S.2)...iT(S n )-v. For example, if one considers dis- 



(8) 



Therefore, if one has sufficient knowledge of the 
underlying noise distribution and the probability 
distributions due to individual sources, then the 
bin counts, rij, can be determined from the ob- 
served pixel distribution with standard nonlinear 
X 2 -minimization techniques applied in the Fourier 
domain. The fitting must also include a parame- 
ter SS for the true underlying flux corresponding 
to map zero values. The differential source counts 
are directly proportional to the bin counts and are 
expressed as dN/dS(Si) — rii/AAi, in terms of the 
fitted map area A and the bin width Aj. 

The pixel distributions arising from sources re- 
lies on the knowledge of the point spread function 
(PSF), which in our case is the closely Gaussian 
shape of the LABOCA beam (Fig. [2]), which has 
been smoothed and large-scale filtered exactly like 
the input map. The underlying noise distribution 
was obtained directly from the jackknifed maps 
(Sect. 13.1) ). which were also smoothed and large- 
scale filtered identically to the input map. 

It is practical to apply the method on the 
signal-to-noise rather than the flux distribution 
of the map, as here the noise distribution is the 
narrowest when coverage is not uniform. This 
helps to limit the unwanted 'smearing' of the ana- 
lyzed probability distribution, thereby improving 
sensitivity for determining the underlying source 
counts. In our case, due to the highly uniform 
coverage in the central 30' x 30' area of our map 
used for determining source counts, the choice be- 
tween flux or S/N distribution fitting is less crit- 
ical. Nevertheless, we analysed the source counts 
on the beam-smoothed S/N image. 

Our choice of a non-linear fitting alg orithm was 
based on the downhill simplex method (jPress et al. 

To use high bin resolution, necessary 
given the expected steepness of the counts (e.g 



Coppin et al. l2006h . and at the same time limit 



the number of fitted parameters to a handful (in 
order to minimize covariances between them and 
to obtain precise fit values for each of these) we 
fitted common source-counts laws rather than the 
individual bin counts n, . 
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Fig. 10. — Left: Differential source counts for the ECDFS. The colored lines show the results from the 
P(D) analysis with functional parameters as given in TableHJ The black data points show the results from 
a Bayesian approach to estimate the source counts fro m the source catalo g. The grey data points are the 
differential number counts from the SHADES survey (|Coppin etal.ll2006f) , the grey line shows their best 
fitting Schechter function. The dashed grey line is the SHADES Schechter function with N' scaled to fit our 
Bayesian source counts. Note that this source count function does not reproduce the observed flux density 
histogram of the map well (Fig. [5]) . Right: Cumulative number counts for the ECDFS compared to other 
studies. The black data points represent the direct sum of the differential counts shown in the left part of 
the figure, the solid lines are i ntegrals over the resu lts of the P(D) analysis, the dashed line shows the best 
fitting Schechter function from lCoppin etalj (|200d ) for SHADES. 



We have fitted four different source count mod- 
els to the observed histogram: a single and a bro- 
ken power law, a Schechter function and a power 
law with constant counts at the faint end of th e 
distribution as suggested by iBarger et al. ( 19991 ). 
All models give comparable fits to the observed 
flux histogram with reduced % 2 values about 1. 
The parameters for each model are given in Ta- 
ble [TJ As only the Barger source count model 
yields finite counts for flux densities approaching 
zero we also give the cutoff flux density derived 
from the fitting together with the implied total 
extragalactic 870/im background light (EBL) con- 
tribution in the Table. Fig. [9] shows an example of 
the simulated flux histogram for the best fitting 
Schechter function in comparison to the observa- 
tions. Within the 3<r Poisson errors there is no 
significant difference between the model and the 
observations. The differential source counts are 
listed in Table [1] and are shown in Fig.fTUl 



3.3.2. Direct estimate of the differential source 
counts 

To compute the differential source counts di- 
rectly from the source catal og, we have u sed a n 



approach similar to that of ICoppin et al. ( 20061 ). 



To determine the posterior flux density distribu- 
tion for each source we have used the power law 
function derived from our P(D) analysis (TableQ]) 
as a prior for the Bayesian deboosting of the ob- 
served flux density. The differential source counts 
were calculated from the deboosted fluxes taking 
the Gaussian likelihood that a source falls into 
a given flux bin as well as the completeness of 
this bin into account. The resulting differential 
source counts are shown together with the best 
fitting models from our P(D) analysis in Fig. [10] 
(left). In the figure we show for comparison the 
differ ential source count s from the SHADES sur- 
vey (jCoppin et al. 20061 ) which is comparable in 
size, but has a noise level ~ 2 higher than LESS. 
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Table 1 

Best fitting parameters of the differential source counts to the observed flux 

histogram. 
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For comparison to previous work, we derived 
the cumulative source counts by directly summing 
over the differential source counts derived above. 
The cumulative source counts are shown in com- 
parison to other studies in Fig. [TO] (right). In this 
figure we also show the integrals over the functions 
fitted by our P(D) analysis (TableHJ. 

3.4. Two point correlation 

We have investigated the clustering prop- 
erties of the SMGs in the ECDFS by means 
of an angular two point correlation function. 
w(6 ) and its uncer t ainty was computed using 
the (Landv fc Szalavl {l993) estimator. The ran- 
dom catalog was generated from the same sim- 
ulations we used for our completeness estimate 
(Sect. l3~2~2|) . To generate random positions of the 
sources we used the LINUX ra ndom number gen- 
" 120061 ). The angular two 



erator (Gutt erman et al 



point correlation is presented in Fig. [IT] We de- 
tect positive clustering for angular scales below 
~ 1', although only the smallest angular scale 
(20" — 50" bin) shows statistically significant clus- 
tering (3.4 a). For comparison to other studies we 
fit the angular correlation by a single power law 
using 

w{9) = A w (^-tO - CO, (9) 

where C accounts for the bias to lower values of 
the observed compared to th e true correlation (see 
e.g. Brainerd fc Smail I ll998h . As our data are too 
noisy to fit all three parameters, we fixed 7 to 1.8 
which has been used in many other studies (e.g 



i Daddi et afll200d iFarrah et al.ll2006t lHartlev et al 
l2008h . This yields A w = 0.011 ± 0.0046 and 
C = 12.4 ± 2.5 or a characteristic clustering an- 
gle of 6 = 14" ± 7". We also calculated C di- 
rectly from our ra ndom catalog using Eq. 22 from 
IScott et all (|2006l ) and assuming 7 = 1.8. This 
yields C=4.5 and A w = 0.007±0.004 (0 O = 7"±5") 
for a single parameter fit of Eq.[9] to our data. 
These numbers demonstrate that the it remains 
difficult to derive the strength of the SMG cluster- 
ing from our data but also that higher spatial res- 
olution would greatly help to identify close SMG 
pairs for a better determination of the clustering. 

4. Discussion 

4.1. Source Counts 

The source counts derived from our P(D) anal- 
ysis show a reasonable agreement with the di- 
rect estimate from the source catalog (Fig.[TOJ. 
We find, however, that the direct counts yield 
significantly lower number counts for flux densi- 
ties of ~ 3 mJy and somewhat higher counts for 
~ 7mJy compared to the P(D) counts. The di- 
rect counts therefore suggest a deviation from a 
single power law with a brea k between 5 — 7mJy , 
similar to the break found bv lCoppin etalJ J2006h 



for th e SHADES survey and by iKnudsen et al 



(|2008D for the Leiden SCUBA Lens Survey. This 
is shown in Fig. [TO] (left) where we plot the best 
fitting Schechter fun ction for the SHADES survey 
( Coppin et al. 20061 ) with N' adjusted to fit our 
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Fig. 11. — Observed angular two point correlation 
function. The solid curve shows a power law fit to 
the data, the dotted line is the corresponding func- 
tion derive d for SMGs select ed from a variety of 
surveys by IScott et al.1 (|2006l ).The dashed vertical 
line shows the size of the LABOCA beam. 



direct counts. This comparison suggests a similar 
shape (but different normalization) of the source 
counts between LESS and SHADES. A compari- 
son of the resulting flux density histogram of this 
model (P(D) like diagram) to the observations, 
however, clearly shows a significant deviation of 
the pixel counts for flux densities below ~ 4mJy 
(Fig.®. 

The most likely explanation for this discrepancy 
is that the direct source count method does not 
take multiple sources in the beam into account, 
but assumes that the deboosted flux density arises 
from a single source. However, from our Monte 
Carlo simulations we find that this assumption is 
poorly justified for LESS. For differential counts 
following a single power law as given in Table[T] 
we find that 25% of the extracted sources in the 
simulations are in fact multiple sources that fall 
too close together (and/or have too poor signal 
to noise) to be identified as double sources. As 
these multiple sources are recovered as a single 
source in the extraction, this naturally leads to an 
underestimate of the faintest sources while over- 
estimating counts at higher flux densities. This 
explains the observed differences between both 
methods and therefore we are confident that the 
P(D) analysis (which does take source overlap into 
account) yields more reliable results on the shape 



of the submm source counts. We note that for 
the source counts derive d from the SHADES sur- 
vey Coppin et al. ( 20061 ) used a correction factor 
determined from Monte-Carlo simulations to take 
the underestimate of the faintest sources into ac- 
count. 

In any case, both methods show consistently 
that submm number counts in the ECDFS are sig- 
nificantly lower for flux densities above ~ 3 mJy 
compared to any other deep-fields observed at 
850 fim so far. Using the SHADES number counts 
for comparison (which are representative also for 
other fields observed with SCUBA, see Fig. fTO]) 
we find that the cumulative source counts of the 
ECDFS are smaller by a factor of ~ 2 for flux 
densities above 3 mJy. This is in line with results 
from optical/NIR surveys which revealed that sev- 
eral restframe-optical populations are underabun- 
dant in the CDFS comp ared to other deep fields: 
van Dokkum et al. ( 20061 ) showed that massive K- 
band selected galaxie s at redshift > 2 are un- 
derabundant by 60%. iMarchesini et al. ( 2007 ) re- 
ported a mild underd ensity of z > 2.5 DRGs and 
Dwellv fc Pagd (|2006[ ) reached similar conclusions 
for optically bright AGNs at high redshifts. These 
studies only covered the central 15' x 15' of the field 
where the submm source density is even lower than 
in the full field (see Sect. [4~3] and Fig.[12]). Simi- 
lar underdensities have also been reported in the 
30' x 30' ECDFS for BzKs dBlanc et al l 120081) as 



well as for X-ray sel ected sources ([Lehmer et al 



20051 : iLuo et al.ll2008l ). The central region of the 



ECDFS is also under-dense in faint (Si.4gh z ~ 
40yuJy) radio sources (E. Ibar priv. comm.). We 
note, however, that the SMC underdensity is only 
observed for bright sources. Due to the steep 
slope and the lack of evidence for a flattening of 
the source counts from our P(D) analysis for low 
fluxes, our number counts become consistent with 
results derived for the faint end of the SMG pop- 
ulation in gravita tionally lensed fields at flux den- 
sities of ~ 1 m Jy (jSmail et alj|2002i iKnudsen et al 



2008). 



4.2. Clustering 

Evidence for strong clustering of SMGs first 
emerged from overdensities of SMGs in close 
proximity to other high-redshift objects (e.g. 
Ivison et a 2000: Chapman et al. 2001; Stevens, et al 



2003. 2004 iKneib et alj 12004 iBeelen et all 120081 ) 
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Direct measurements of 3 an d 2-D clustering o f 



SMGs have been p resented by Blain et al 



Greve et all (|2004h and IScott et all (|2009 1. Due 



(2004), 



to the small number of sources typically involved 
in these studies, the significance of the cluster- 
ing amplitude, however, remains marginal. Al- 
though the number of source detected by us in 
the ECDFS only gives a small improvement in 
terms of signal to noise compared to the previous 
measurements, our detection of strong clustering 
in an independent submm survey greatly improves 
the reliability of SMG clustering results. 

Our characteristic clustering angle of 9q ~ 
5 — 15" is smaller than the a n gle o f 9o ~ 
40 - 50" derived by IScott et all (|2006l ) but in 



agr eement with the predict ions of merger mod- 
els ( van Kampen et al. 2005h . which predict clus- 
tering scales between 5" (for a hydrodynamical 
model) and 20" (for a high mass merger model). 
The differe n ce be tween our results and those by 
IScott et all (|2006l ) can most likely be explained 
by the small significance of the clustering signal 
in both studies. From Fig.[TT]it can be seen that 



the powerlaw fit by IScott et al. ( 2006f) 



is consis- 
tent with our most significant data point at an 
angular distance of 40" and that our smaller clus- 
tering angle mainly results from the small angular 
correlation for distances b e tween 1' — 2' in our 
data. However I Greve et"aL ( 2004h also presented 
evidence for an excess of SMG pairs with a typical 
separation of 23" based on 1.2 mm observations 
which appears to support our findings. We note 
that the angular correlation function should not 
be affected by the LSS filtering of our map because 
it refers to the angular seperation between point 
sources which remains unaffected by the filtering. 

A comparison of the clustering of SMGs to 
other high rcdshift source populations has the 
potential to shed more light on the evolution of 
SMGs and to investigate to what kind of sources 
SMGs evolve once the gas has been consumed. 
The clustering, h owever, is expecte d to evolve 
with redshift (e.g. Farrah et al. 2006f ) and a con- 
clusive comparison of the clustering relies on a 3- 
D clustering analysis which requires knowledge of 
the redshift distribution of our SMGs. Assum- 
ing the same redshift dist ribution for our SM Gs 
as the distribution used bv lFarrah et al. I (l2006h (a 
Gaussian centered at z=2.5 with a FWHM of 1.2, 
similar to the redshift distribution of SMGs de- 



rived bv lChapman et al" ( 20051 )) our data suggests 
a correlation length of ro = 10 ± 6/i _1 Mpc and 
7'o = 13 ±6 h~ x Mpc for our single and two param- 
eter fit, respectively. These values are larger than, 
but consistent within the errors with, the correla- 
tion length oforn = 6.9 ± 2.1 h' 1 Mpc derived by 



Blain et al 



2004h for SMGs and in good agree- 
ment with the correlation length of ro = 14.4 ± 
2. hT 1 Mp c derived for 24 ^m selected ULIRGs at 



2-3 jFarrah et al.ll2006f) 



Most models for the evolution of overdensitics 
over time predict an increasing correlation length 
for decreasing re dshift (for a collection of evolution 



models see e.g. lOverzier et al.l 120031 ) . This sug 



gests that the successors of SMGs could be asso- 
ciated with clusters of galaxies (ro ~ 20 ft. -1 Mpc) 
at the present epoch. Comparing our clustering 
strength to the correlation length of dark mat- 
ter (DM) halos as a function of redshift (e.g. 
iMatarrese et al. 1997 ) furthermore suggests that 
SMGs reside in more massive (~ 10 13-14 M Q ) DM 
halos than other high-z source populations such 
as LBGs and QSOs. Given the uncertainties of 
our measurement, the unknown redshift distribu- 
tion of our SMGs and the model dependence of 
the DM clustering these numbers are quite uncer- 
tain but_sjap_p_OTt_the conclusions of prev ious stud- 
ies (|Blain et al.ll2004l : iFarrah et al.ll2006h . 



4.3. Spatial variations of the source counts 

Motivated by the strong clustering detected in 
the distribution of the submm sources we also 
investigated the spatial variations of the source 
counts across the map. To distinguish between 
regions with potential over and underdensities of 
submm sources we have used the integrated ex- 
tragalactic 870^m flux detected in the map as ref- 
erence. The variations were calculated at spatial 
resolution of 5' by convolving the flux map by a 
Gaussian kernel. The resulting flux distribution is 
shown in Fig. [12] in units of Jydeg -2 and reflects 
mainly the large scale distribution of the individ- 
ually detected sources as faint extended emission 
has been removed by our optimal point-source fil- 
tering operation (see Sect. [2]). We have used the 
4.1 Jydeg -2 contour to divide our source catalog 
into two sub-catalogs which we call for convenience 
the 'sparse' and the 'dense' region in the following. 
The contour was chosen to give roughly an equal 
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Jy deg 




Fig. 12. — 870 fim flux density distribution 
smoothed to 5' spatial resolution. The flux den- 
sity is show in Jydeg -2 . The black square shows 
the central 30' x 30' of the ECDFS. The white cir- 
cles mark the location of the 126 sources listed in 
Table [2j Black contours are shown from 1 to 7 Jy 
deg -2 spaced by 1 Jy deg -2 . The white contour is 
the 4.1 Jydeg -2 level that has been used to define 
the SMG over and underdensity region. 



number of detected sources for both catalogs. The 
source selection based on the large scale flux dis- 
tribution emphasizes the difference of the source 
counts between over and underdensities and is no 
longer representative for the counts in a blank field 
survey. 

The source counts were determined in the same 
manner as described in Sect. 13. 31 using a P(D) 
analysis and a direct estimate of the source counts. 
In the top panel of Fig.[T2] we show the observed 
flux histograms for both regions in comparison to 
the fits from the P(D) analysis. For the sparse 
region we find that the source counts are well 
described by a single power law. The best fit- 
ting model yields a normalization at 5 mJy half of 
that for the full field (N 1 = 47 ± 8 deg -2 mJy -1 ) 
and a slope of a — 3.6 ± 0.3 - steeper than the 
counts for the full field {a = 3.2 ± 0.2). The 
P(D) analysis suggests that the source counts of 
the dense region are much shallower. For a single 
power law we find good matching parameters of 
N' = 250 ± 20deg -2 mJy -1 and a = 2.9 ± 0.2. 
The fitting, however, suggests that the counts are 



slightly better described by a broken power law 
or a Schechter function. For a broken power law 
the break occurs at ~ 8 mJy with a normaliza- 
tion of N' = 100 ± 8 deg -2 mJy -1 . For the slopes 
we find a = 2.4 ± 0.15 and /? = 4.7 ± 0.6. In 
the bottom part of Fig. [13] we compare the inte- 
grated source counts determined from the direct 
source count estimate and the P(D) analysis for 
both regions to the counts of other 870/im surveys. 
Fig. [TBI shows that the source counts for the dense 
region are in good agreement with results of the 
SHADES and other submm surveys. This com- 
parison demonstrates that 1) the surface density 
of submm sources is not constant, but changes by 
a factor ~ 3 on angular scales of ~ 10' and 2) the 
sampling variance is not due to a simple scaling 
of the number counts but associated with changes 
of the shape of the source counts. The latter find- 
ing suggests that the sampling variance is not due 
to (weak) lensing by foreground mass distributions 
but intrinsic to the distribution of submm sources. 

4.4. Contribution to the EBL 

The P(D) analysis of the differential source 
counts also provides an estimate on the integrated 
870 /im background light that can be compared to 
the interpolated EBL at subm m wavelength from 
COBE FIR AS of ~ 45 Jydeg -2 (jPuget et alJl994 
19981) . Depending on the assumed 



iFixsen et al 



underlying source distribution we recover an EBL 
contribution of 29 - 32 Jydeg -2 for the ECDFS 
(see Table[TJ). We note that this range is a lower 
limit to the true underlying 870 flux as the con- 
tribution of very faint sources, which are expected 
to form an almost uniform distribution across the 
field at the spatial resolution of LABOCA, will be 
removed from our data in the correlated noise re- 
moval steps. In comparison to the values measured 
by COBE we therefore detect > 65 - 70% of the 
EBL for sources brighter than ~ 0.5 mJy (which 
corresponds to a typical lower flux cutoff in our 
P{D) analysis). These numbers demonstrate that 
it is unlikely that the integrated EBL level in the 
ECDFS is significantly lower than in other parts 
of the sky. This implies that the observed factor 
2 underdensity of submm source relative to other 
deep fields is restricted to ULIRGs with far in- 
frared luminosities of > 2 x 10 12 L Q (assuming no 
lensing and z > 0.5 for the bulk of our sources), 
while more typical star forming galaxies, which 
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S„ [mJy] 

Fig. 13. — Variation of the source counts with 
source surface density: Top: The black and the 
red histograms show the pixel flux distribution for 
the dense and sparse region fFig. [T2")) . The green 
curves show the results from our P(D) analysis. 
All histograms are normalized to a peak value of 
1. Bottom: Integrated source counts in both sub- 
regions. The data points are the source counts 
derived from the direct estimate, the black lines 
the results from the P(D) analysis. The grey line 
shows the SHADES source counts, the grey area 
the source counts derived by other 870 fim studies 
(see FigfTOl right). 

dominate the EBL, are not underabundant in the 
field. This conclusion is also supported by the in- 
tensity of the EBL we find for the two subregions 
discussed in Section FOl where our P{D) analysis 
does not yield a significant difference between the 
two fields. 

5. Summary and Conclusions 

We have presented a deep 870 survey of the 
ECDFS using LABOCA on the APEX telescope at 



Llano de Chajnantor in Chile. This is the largest 
contiguous deep submm survey to date. Our map 
has a highly uniform noise level across the full 30' x 
30' field of 1.2 mJy beam -1 and our survey is > 
95% complete for sources down to a flux limit of 
6.5 mJy. Our main findings are summarized as 
follows: 

• At the (beam smoothed) spatial resolution 
of 27" of our survey we find that the map's 
noise level is affected by confusion noise 
arising from faint, individually undetected 
SMGs. From the rms noise as a function of 
integration time we derive a confusion noise 
of a c w 0.9 mJy beam -1 . 

• We identify 126 submm sources in a search 
area of 1260arcmin 2 above a signal to noise 
threshold of 3.7 cr, which corresponds to an 
expected false detection rate of 5 sources. 

• We have determined the differential and in- 
tegrated source counts using a P(D) analy- 
sis and an estimate based on our source cat- 
alog. Both results are in reasonable agree- 
ment and show that SMGs in the ECDFS 
are underabundant by a factor of ~ 2 for 
sources brighter than 3 mJy compared to the 
average of previous surveys. Under the as- 
sumption that the bulk of the sources are 
at z > 0.5, this implies an underdensity of 
ULIRGs with Lfir > 2 x 10 12 L q compared 
to other blank fields that have been observed 
in the submm. The source counts are well 
described by a single power law with a slope 
of a = 3.2 ±0.2. 

• We derive the angular two-point correlation 
function for the SMGs and find clustering on 
angular scales < 1' with a significance up to 
3.4 cr. Assuming a power law dependence for 
the correlation function we derive a cluster- 
ing amplitude of A w = 0.011 ± 0.0046 or a 
characteristic angular scale of 9q = 14" ± 7" 
for 7 = 1.8. Assuming a redshift distribu- 
tion similar to that observed for spectroscop- 
ically confirmed SMGs, we derive a correla- 
tion length of ro = 13 ± 6 hr 1 Mpc, some- 
what larger than previous estimates of the 
3-D clustering of SMGs but in agreement 
with the clustering derived for 24 (im se- 
lected ULIRGs. 
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• We have investigated for the first time the 
spatial variations of the SMG source counts. 
We find that the differential source counts in 
regions with an ovcrdcnsity of SMGs have 
a different shape compared to those with 
underdensities. While the counts in un- 
derdense regions are well fitted by a single 
power law with a slope of a = 3.6 ± 0.3, 
the counts in the overdensities are signifi- 
cantly shallower with a — 2.9 ± 0.2. The 
counts in the overdensities are slightly bet- 
ter described by a broken power law or a 
Schechter function. For flux densities below 
8mJy we find a slope of a = 2.4 ± 0.f5, for 
sources above this limit the counts are much 
steeper with (3 = A.l ± 0.6. This may indi- 
cate an intrinsic turn-over in the underlying 
luminosity function placing an upper limit 
on the FIR luminosity. 

• The integrated 870 /xm flux density derived 
from our survey is > 29 — 32 Jydeg -2 for 
sources brighter than ~ 0.5 mJy which cor- 
responds to > 65 — 70% of the extragalac- 
tic background light estimated from COBE 
measurements. We do not find a signifi- 
cant difference of the quantity between SMG 
over- and underdensities. We conclude that 
ECDFS is undcrabundant of ULlRGs but 
not of more typical star forming systems 
with lower FIR luminosities, which dominate 
the extragalactic background light. 



We would like to thank the APEX staff for 
their aid in carrying out the observations. APEX 
is operated by the Max-Planck-Institut fur Ra- 
dioastronomie, the European Southern Observa- 
tory, and the Onsala Space Observatory. IRS, 
KEKC and RJI acknowledge support from STFC. 
JSD acknowledges the support of the Royal Soci- 
ety and of STFC. 



20 



Table 2 

870/iM LABOCA source catalog of the ECDFS. 
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.8 


9.1 ± 1.7 


8.1 ± 


1 


.8 


5 


.4 


0.0 


2 


LESSJ033149.2-280208 


03 


31 


49.15 


-28 


02 


08 


.7 


6.9 ± 1.3 


6.4 ± 


1 


.4 


5 


.4 


0.0 


2 


LESS J033336.0-275347 


03 


33 


36.04 


-27 


53 


47 


.6 


7.3 ± 1.4 


6.7 ± 


1 


.5 


5. 


.3 


0.0 





LESSJ033310.2-275641 


03 


33 


10.20 


-27 


56 


41 


.5 


6.4 ± 1.2 


6.0 ± 


1 


.3 


5. 


.2 


0.0 





LESSJ033144.9-273435 


03 


31 


44.90 


-27 


34 


35 


.4 


6.8 ± 1.3 


6.2 ± 


1 


.4 


5. 


.2 


0.0 





LESSJ033246.7-275120 


03 


32 


46.74 


-27 


51 


20 


.9 


6.4 ± 1.2 


5.9 ± 


1 


.3 


5. 


.2 


0.0 


2 


LESSJ033110.5-275233 


03 


31 


10.47 


-27 


52 


33 


.2 


8.7 ± 1.7 


7.6 ± 


1 


.9 


5 


.2 


0.0 


1 


LESSJ033231.0-275858 


03 


32 


31.02 


-27 


58 


58 


.1 


6.4 ± 1.2 


5.8 ± 


1 


.4 


5. 


.1 


0.0 


2 


LESSJ033307.0-274801 


03 


33 


07.00 


-27 


48 


01 


.0 


6.4 ± 1.3 


5.9 ± 


1 


.4 


5. 


.1 


0.0 





LESSJ033131.0-273238 


03 


31 


30.96 


-27 


32 


38 


.5 


7.5 ± 1.5 


6.7 ± 


1 


.6 


5 


.1 


0.0 





LESSJ033225.7-275228 


03 


32 


25.71 


-27 


52 


28 


.5 


6.3 ± 1.2 


5.8 ± 


1 


.4 


5 


.1 


0.0 


1 


LESS J033336.8-273247 


03 


33 


36.80 


-27 


32 


47 


.0 


8.2 ± 1.6 


7.2 ± 


1 


.8 


5 


.1 


0.0 
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T TTIOC Tnoorir / ' 

LEbb J033256 


.0 


-273317 


03 


32 


56 


.00 


-27 


33 


17 


.7 


7. 


.0 


± 


1 


.4 


6. 


.3 


± 


1 


.5 


5 


.1 





.0 





LEbb J 033237 


.8- 


-273202 


03 


32 


37 


.77 


-27 


32 


02. 


.0 


7. 


.7 


± 


1 


.5 


6 


.8 


± 


1 


.7 


5 


.1 





.0 





LEbb J033124 


.5- 


-275040 


03 


31 


24 


.45 


-27 


50 


40 


.9 


6. 


6 


± 


1 


.3 


5. 


.9 


± 


1 


.4 


5 


.1 





.0 


1,2 


LEbb J033141 


.2- 


-274441 


03 


31 


41 


.15 


-27 


44 


41 


.5 


6. 


.1 


± 


1 


.2 


5. 


.6 


± 


1 


.3 


5 


.0 





.0 


1,2 


T TTIfP TAOOI A A 

LEbb J 033144 


.8- 


-274425 


03 


31 


44 


.81 


-27 


44 


25 


.1 


6. 


.2 


± 


1 


.2 


5. 


.6 


± 


1 


.3 


5 








.0 


2 


T TTIOC TAOH AO 

LEbb J033I28 


.5- 


-275601 


03 


31 


28 


.51 


-27 


56 


01 


.3 


6. 


.2 


± 


1 


.3 


5. 


.6 


± 


1 


.4 


4. 


.9 





.0 


1 


LEbb J033159 


.1- 


-275435 


03 


31 


59 


.12 


-27 


54 


35 


.5 


6. 


.2 


± 


1 


.3 


5. 


.6 


± 


1 


.4 


4. 


.9 





.0 


2 


LEbb J 033243 


.6 


-273353 


03 


32 


43 


.61 


-27 


33 


53 


.6 


6 


.8 


± 


1 


.4 


6. 


.0 


± 


1 


.5 


4, 


.9 





.0 


2 


LEbb J033302 


.2- 


-274033 


03 


33 


02 


.20 


-27 


40 


33 


.6 


6. 


.1 


± 


1 


.2 


5. 


.5 


± 


1 


.4 


4. 


.9 





.0 





LEbb J033I53 


.2 


-273936 


03 


31 


53 


.17 


-27 


39 


36 


.1 


6. 


.0 


± 


1 


.2 


5 


.4 


± 


1 


.4 


4. 


9 





.0 


1 


LEbb J033I52 


.0-275329 


03 


31 


51 


.97 


-27 


53 


29 


.7 


6. 


.1 


± 


1 


.3 


5 


.5 


± 


1 


.4 


4. 


9 





.0 





LEbS J033225 


.8- 


-273306 


03 


32 


25 


.79 


-27 


33 


06 


.7 


6 


7 


± 


1 


.4 


5 


.9 


± 


1 


.6 


4. 


.8 





.0 





LEbb J033303 


.9 


-274412 


03 


33 


03 


.87 


-27 


44 


12 


.2 


6 





± 


1 


.3 


5. 


.3 


± 


1 


.4 


4. 


.8 





,0 





LEbb J033317 


.5- 


-275121 


03 


33 


17 


.47 


-27 


51 


21 


.5 


5. 


.8 


± 


1 


.2 


5. 


.2 


± 


1 


.4 


4. 


.8 





.1 


1 


LEbb J033245 


.6 


-280025 


03 


32 


45 


.63 


-28 


00 


25 


.3 


5 


9 


± 


1 


.2 


5. 


.2 


± 


1 


.4 


4. 


.7 





.1 





LEbb J033236 


.4 


-273452 


03 


32 


36 


.41 


-27 


34 


52 


.5 


6. 


.1 


± 


1 


.3 


5. 


.4 


± 


1 


.5 


4. 


.7 





.1 


2 


LEbb J033308 


.5- 


-280044 


03 


33 


08 


.46 


-28 


00 


44 


.3 


6. 


.0 


± 


1 


.3 


5 


.3 


± 


1 


.4 


4. 


7 





.1 


2 


T TOO TflOOOfH 


.0- 


-280025 


03 


32 


01 


.00 


-28 


00 


25 


.6 


5. 


.8 


± 


1 


.2 


5 


.1 


± 


1 


.4 


4. 


.7 





.1 


1 


T TTIfP TAOOAPA 

LEbb J033252 


.4 


-273527 


03 


32 


52 


.40 


-27 


35 


27 


.7 


5. 


9 


± 


1 


.3 


5. 


.2 


± 


1 


.4 


4. 


.7 





.1 





T 1 ^ I 1 i ' TAOOOOI 

LEbb J033331 


.7 


-275406 


03 


33 


31 


.69 


-27 


54 


06 


.1 


6 


1 


± 


1 


.3 


5. 


.3 


± 


1 


.5 


4. 


.7 





.1 





T TOO TAOOO/IO 

Li^bo 31)63243 


.3 


-275517 


03 


32 


43 


.28 


-27 


55 


1 1 


.9 


5. 


9 


i 

± 


1 


.3 


5 


.2 


± 


1 


A 

.4 


A 

4 


-7 
. 1 





.1 





T TTIfP TAOOAOO 

LEbb J033233 


.4- 


-273918 


03 


32 


33 


.44 


-27 


39 


18 


.5 


5. 


.8 


± 


1 


.3 


5. 


.1 


± 


1 


.4 


4. 


7 





.1 


1 


T TTiO O TAOOI O A 

LEbb J033134 


.3 


-275934 


03 


31 


34 


.26 


-27 


59 


34 


.3 


5. 


.7 


± 


1 


.2 


5 


.0 


± 


1 


.3 


4. 


.7 





.1 


1 


T TTiOC TAOOI A A 

LEbb J033144 


.0- 


-273832 


03 


31 


43 


.97 


-27 


38 


32 


.5 


5. 


.7 


± 


1 


.2 


5. 


.0 


± 


1 


.4 


4. 


.6 





.1 


2 


T TAOOOA/"" 

LEbb J033306 


.3 


-273327 


03 


33 


06 


.29 


-27 


33 


27 


.7 


6 


6 


± 


1 


.4 


5. 


.6 


± 


1 


.6 


4. 


.6 





.1 





T TTiOO TAOOA/1A 

LEbb J 033240 


.4- 


-273802 


03 


32 


40 


.40 


-27 


38 


02. 


.5 


5. 


.7 


± 


1 


.2 


5 


.0 


± 


1 


.4 


4. 


.6 





.1 


1 


LEbb J033229 


.3- 


-275619 


03 


32 


29 


.33 


-27 


56 


19 


.3 


5. 


8 


± 


1 


.3 


5 


.1 


± 


1 


.4 


4. 


.6 





.1 





T TTiOO T AO O O AA 

LEbb J033309 


.3 


-274809 


03 


33 


09 


.34 


-27 


48 


09 


.9 


5. 


.8 


± 


1 


.3 


5. 


.1 


± 


1 


.4 


4. 


.6 





.1 





LEbb J033126 


.8- 


-275554 


03 


31 


26 


.83 


-27 


55 


54 


.6 


5. 


.8 


± 


1 


.3 


5. 


.1 


± 


1 


.4 


4. 


.6 





.1 


2 


T TT^OO TAOOOOA 

LEbb J033332 


.7-275957 


03 


33 


32 


.67 


-27 


59 


57 


.2 


6. 





± 


1 


.3 


5 


.1 


± 


1 


.5 


4. 


.5 





1 





T TTiOC TAOOI P - T 

LEbb J033157 


.2- 


-275633 


03 


31 


57. 


.23 


-27 


56 


33 


.2 


5. 


.6 


± 


1 


.3 


4. 


.8 


± 


1 


.4 


4. 


.4 





.2 





T TTIOC TAOOO/1A 

LEbb J 033340 


.3- 


-273956 


03 


33 


40 


.30 


-27 


39 


56 


.9 


6 


2 


± 


1 


.4 


5 


.1 


± 


1 


.7 


4. 


.4 





.3 


1 


LEbb J 033221 


.3 


-275623 


03 


32 


21 


.25 


-27 


56 


23 


.5 


5. 


.5 


± 


1 


.3 


4. 


.7 


± 


1 


.4 


4. 


.4 





.3 


2 


T TTiOC TAOOI A A 

LEbb J033142 


.2 


-274834 


03 


31 


42 


.23 


-27 


48 


34 


.4 


5. 


.4 


± 


1 


.2 


4. 


.6 


± 


1 


.4 


4. 


,4 





.3 


2 


LEbb J033127. 


.5- 


-274440 


03 


31 


27 


.45 


-27 


44 


40 


.4 


5. 


.7 


± 


1 


.3 


4. 


.8 


± 


1 


.5 


4. 


.4 





.3 





LEbSJ033253. 


.8- 


-273810 


03 


32 


53 


.77 


-27 


38 


10 


.9 


5. 


.4 


± 


1 


.2 


4. 


.5 


± 


1 


.4 


4. 


.4 





.3 





LEbbJ033308 


.9 


-280522 


03 


33 


08 


.92 


-28 


05 


22 


.0 


6. 


.7 


± 


1 


.5 


5. 


.3 


± 


1 


.8 


4. 


.4 





.3 





LEbbJ033154 


.2- 


-275109 


03 


31 


54 


.22 


-27 


51 


09 


.8 


5. 


.5 


± 


1 


.3 


4. 


.6 


± 


1 


.4 


4. 


3 





.4 





LEbbJ033110 


.3 


-274503 


03 


31 


10 


.28 


-27 


45 


03 


.1 


8, 


2 


± 


1 


.6 


6 


.0 


± 


2 


.4 


4. 


3 





.4 


1 


LEbSJ033114 


.9 


-274844 


03 


31 


14 


.90 


-27 


48 


44 


.3 


6 


.5 


± 


1 


.5 


5. 


.1 


± 


1 


.8 


4. 


3 





.4 





LEbbJ033251 


.1- 


-273143 


03 


32 


51 


.09 


-27 


31 


43 


.0 


6. 


.7 


± 


1 


.6 


5 


.3 


± 


1 


.9 


4. 


3 





.5 





LEbbJ033155 


.2 


-275345 


03 


31 


55 


.19 


-27 


53 


45 


.3 


5. 


.4 


± 


1 


.3 


4. 


.5 


± 


1 


.4 


4, 


3 





.5 





LEbSJ033248. 


.4 


-280023 


03 


32 


48 


.44 


-28 


00 


23 


.8 


5. 


.3 


± 


1 


.2 


4. 


.4 


± 


1 


.4 


4. 


.3 





.5 





LEbSJ033243 


.7-273554 


03 


32 


43 


.65 


-27 


35 


54. 


.1 


5 


.4 


± 


1 


.3 


4. 


.5 


± 


1 


.5 


4, 


.2 





.6 





LEbSJ033135 


.3 


-274033 


03 


31 


35 


.25 


-27 


40 


33 


.7 


5. 


.3 


± 


1 


.3 


4, 


.4 


± 


1 


.4 


4. 


2 





.6 





LEbSJ033138 


.4 


-274336 


03 


31 


38 


.36 


-27 


43 


36 


.0 


5. 


.2 


± 


1 


.2 


4. 


.3 


± 


1 


.4 


4. 


2 





.6 





LEbbJ033110 


.8- 


-275607 


03 


31 


10 


.84 


-27 


56 


07 


.2 


6. 


9 


± 


1 


.7 


5. 


.2 


± 


2 


.0 


4. 


2 





.6 





LEbb J033307 


.3 


-275805 


03 


33 


07 


.27 


-27 


58 


05 


.0 


5. 


.3 


± 


1 


.3 


4. 


.4 


± 


1 


.4 


4, 


2 





.7 
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LESS J033241 


.7 


-275846 


03 


32 


41 


.74 


-27 


58 


46 


.1 


5. 


.2 


± 


1 


.3 


4. 


.3 


± 


1 


.4 


4. 


.2 





.7 





LESSJ033313 


.0- 


-275556 


03 


33 


13 


.03 


-27 


55 


56 


.8 


5. 


.2 


± 


1 


.2 


4. 


.3 


± 


1 


.4 


4 


.2 





.7 





LESSJ033313 


.7 


-273803 


03 


33 


13 


.65 


-27 


38 


03 


.4 


5. 


.1 


± 


1 


.2 


4. 


.2 


± 


1 


.4 


4 


.2 





.8 





LESSJ033130 


.2 


-275726 


03 


31 


30 


.22 


-27 


57 


26. 


.0 


5. 


.1 


± 


1 


.3 


4 


.2 


± 


1 


.4 


4 


.1 


1 


.0 





LESSJ033251 


.5- 


-275536 


03 


32 


51 


.45 


-27 


55 


36 


.0 


5. 


.3 


± 


1 


.3 


4 


.3 


± 


1 


.4 


4 


.1 


1 


.0 





LESSJ033111 


.3- 


-280006 


03 


31 


11 


.32 


-28 


00 


06. 


.2 


6. 


.4 


± 


1 


.6 


4 


.8 


± 


1 


.9 


4. 


.1 


1 


.0 





LESSJ033151. 


.5- 


-274552 


03 


31 


51 


.47 


-27 


45 


52 


.1 


5. 


.1 


± 


1 


.3 


4 


.2 


± 


1 


.4 


4. 


1 


1 


.1 


1 


LESSJ033335 


.6 


-274020 


03 


33 


35 


.61 


-27 


40 


20 


.1 


5. 


.4 


± 


1 


.3 


4. 


.3 


± 


1 


.5 


4. 


.1 


1 


.2 





LESSJ033325 


.4 


-273400 


03 


33 


25 


.35 


-27 


34 


00 


.4 


5. 


.5 


± 


1 


.3 


4 


.3 


± 


1 


.5 


4. 


.1 


1 


.2 





LESSJ033258 


.5- 


-273803 


03 


32 


58 


.46 


-27 


38 


03 


.0 


4. 


.9 


± 


1 


.2 


4 


.0 


± 


1 


.4 


4. 


.1 


1 


.3 





LESSJ033115 


.8- 


-275313 


03 


31 


15 


.78 


-27 


53 


13 


.1 


6. 


.0 


± 


1 


.5 


4 


.6 


± 


1 


.7 


4. 


.1 


1 


.3 





LESSJ033140 


.1- 


-275631 


03 


31 


40 


.09 


-27 


56 


31 


.4 


4. 


.9 


± 


1 


.2 


4 


.0 


± 


1 


.4 


4. 


.0 


1 


.4 


2 


LESSJ033130 


.9 


-275150 


03 


31 


30 


.85 


-27 


51 


50 


.9 


5. 


.0 


± 


1 


.3 


4 


.0 


± 


1 


.4 


4. 


.0 


1 


.4 





LESSJ033316 


.4 


-275033 


03 


33 


16 


.42 


-27 


50 


33 


.1 


5. 


.0 


± 


1 


.2 


4 


.0 


± 


1 


.4 


4. 


.0 


1 


.5 





t ttqq thqqqoq 
-L-CjOO JUooozo 


, 1 


-Z 1 410 1 


Do 


•it) 


zo 




-z ( 


1 1 


O i . 


U 


0. 


U 


1 

± 


1 
1 


Q 

.0 


A 
4 


n 
U 


_j_ 
± 


1 
1 


.4 


A 

4. 


n 
U 


1 

1 . 


n 



U 


LESS J033122. 


.6 


-275417 


03 


31 


22 


.64 


-27 


54 


17 


.2 


5. 


.3 


± 


1 


.3 


4 


.1 


± 


1 


.5 


4. 


.0 


1 


.8 





LESSJ033325 


.6 


-273423 


03 


33 


25 


.58 


-27 


34 


23 


.0 


5. 


.2 


± 


1 


.3 


4 


.1 


± 


1 


.5 


4. 


.0 


1 


.9 





LESSJ033249. 


.3- 


-273112 


03 


32 


49 


.28 


-27 


31 


12 


.3 


6. 


.5 


± 


1 


.7 


4 


.6 


± 


2 


.0 


4. 


.0 


2 


.0 





LESSJ033236 


.4- 


-275845 


03 


32 


36 


.42 


-27 


58 


45 


.9 


5. 


.0 


± 


1 


.3 


3. 


.9 


± 


1 


.4 


3. 


.9 


2 


.0 





LESSJ033150. 


.8- 


-274438 


03 


31 


50 


.81 


-27 


44 


38 


.5 


4. 


.9 


± 


1 


.3 


3 


.9 


± 


1 


.4 


3. 


.9 


2 


.4 





LESS J033349 


.7 


-274239 


03 


33 


49 


.71 


-27 


42 


39 


.2 


7. 


.4 


± 


1 


.6 


4. 


.6 


± 


2 


.4 


3. 


.9 


2 


.9 





LESS J033154. 


.4 


-274525 


03 


31 


54 


.42 


-27 


45 


25 


.5 


4. 


.9 


± 


1 


.3 


3. 


.8 


± 


1 


.4 


3. 


.8 


3 


.0 





LESSJ033128. 


.0 


-273925 


03 


31 


28 


.02 


-27 


39 


25 


.2 


5. 


.0 


± 


1 


.3 


3. 


.8 


± 


1 


.4 


3. 


.8 


3 


.1 





LESSJ033121. 


.8- 


-274936 


03 


31 


21 


.81 


-27 


49 


36 


.8 


5. 


.2 


± 


1 


.4 


3. 


.8 


± 


1 


.5 


3. 


.8 


3 


.7 





LESSJ033256. 


.5- 


-280319 


03 


32 


56 


.51 


-28 


03 


19 


.1 


5. 


.1 


± 


1 


.4 


3. 


.8 


± 


1 


.5 


3. 


.8 


3 


.8 





LESSJ033328 


.5- 


-275655 


03 


33 


28 


.45 


-27 


56 


55. 


.9 


4. 


.9 


± 


1 


.3 


3. 


.7 


± 


1 


.5 


3. 


.8 


3 


.8 





LESSJ033333 


.3- 


-273449 


03 


33 


33 


.32 


-27 


34 


49 


.3 


5. 


.2 


± 


1 


.4 


3. 


.8 


± 


1 


.6 


3. 


.8 


3 


.9 





LESSJ033139 


.6 


-274120 


03 


31 


39 


.62 


-27 


41 


20. 


.4 


4. 


.7 


± 


1 


.2 


3. 


.6 


± 


1 


.5 


3. 


.8 


4. 


.0 





LESSJ033330 


.9 


-275349 


03 


33 


30 


.88 


-27 


53 


49 


.3 


4. 


.9 


± 


1 


.3 


3. 


.7 


± 


1 


.6 


3. 


.8 


4. 


.2 





LESSJ033203 


.6 


-273605 


03 


32 


03 


.59 


-27 


36 


05 


.0 


4. 


.6 


± 


1 


.2 


3. 


.5 


± 


1 


.4 


3. 


.7 


4. 


.7 





LESSJ033146. 


.0- 


-274621 


03 


31 


46 


.02 


-27 


46 


21 


.2 


4. 


.7 


± 


1 


.3 


3. 


.6 


± 


1 


.4 


3. 


.7 


4 


.7 





LESSJ033209 


.8- 


-274102 


03 


32 


09 


.76 


-27 


41 


02 


.0 


4. 


.7 


± 


1 


.3 


3 


.6 


± 


1 


.4 


3. 


.7 


4 


.9 






a Deboosted fluxes depend on the source count model (see Sect. l3.2!3j) and are only correct in a statistical sense. 

b Expected number of false detections for all sources including the corresponding entry in the Table. 

c The "Remark" entry indicates if a source is detected in the two submaps calculated by splitting the data into half 
(see 13.2.51 . 1,2: the source is detection in both submaps, 1 (2): the source is detected in submap 1 (2), the source 
is not detected in any submap. 
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